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are  explored:   Sinusoidal,  6-function  and  spatially  decaying 
wall  perturbations.   Strong  effects  (resonances)  appear  in 
the  case  of  the  sinusoidal  wall.   In  the  other  two  cases, 
the  resonances  are  not  present,  but  there  are  still  traveling 
disturbances  whose  effects  can  be  important  at  great 
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I.   INTRODUCTION 

The  theory  of  acoustic  wave  propagation  in  bounded 
media  (waveguides  and  waveducts)  has  been  treated  with  a 
certain  intensity  in  the  last  twenty-five  years.   It  is 
possible  to  classify  the  problems  arising  in  these  studies 
in  two  categories.   One  category  involves  the  effects  on 
acoustic  propagation  of  small  inhomogeneities  and  the  other, 
the  effects  on  the  acoustic  field  of  irregularities  on  the 
boundaries  of  the  waveguide  or  waveduct. 

These  two  factors  are  the  most  important  concerning  the 
acoustic  v/ave  propagation  and  M.A.  Isakovich  [4]  in  1956 
said  that  "such  problems,  besides  their  theoretical  interest 
are  also  important  from  a  practical  point  of  view;  thus,  for 
example,  in  super-long-range  waveduct  propagation  of  sound 
in  a  laminarly-inhomogeneous  medium  the  effect  of  even  small 
inhomogeneities  or  corrugation  of  the  boundaries  can  produce 
a  very  great  effect  upon  the  range  of  propagation" . 
M.A.  Isakovich  studied  propagation  in  a  waveduct  with  small 
inhomogeneities  and  in  a  waveduct  with  one  corrugated  wall, 
using  the  method  of  small  perturbations.   This  was  a  natural 
extension  of  the  problem  of  scattering  of  sound  from  one 
single  rough  wall. 

I.e.  Samuels  [5],  in  1958,  studied  the  problem  of 
propagation  of  acoustic  and  electromagnetic  v/aves  in  ducts 
with  two  corrugated  rigid  walls.   He  found,  in  the  case  of 


harmonic  roughness,  "that  waves  are  not  scattered  upstream 
if  the  frequency  of  the  transmitted  signal  exceeds  certain 
critical  values  related  to  the  natural  modes  of  the  waveguide 
and  the  wavelengths  of  the  wall  roughness." 

C.S.  Clay  [6],  in  1963  presented  a  paper  where  the 
effect  of  a  slightly  irregular  boundary  on  waveguide 
propagation  was  studied.   The  results  were  checked  by  an 
experiment  in  shallow  water  and  he  found  that  the  effects  of 
small  roughness  on  the  surface  are  to  increase  the  attenua- 
tion as  a  function  of  range  and  decrease  the  coherence  of 
acoustical  propagation. 

In  more  recent  years  (1971) ,  R.H.  Ebert  [7]  calculated 
the  fluctuation  of  sound  due  to  influence  of  a  standing 
gravity  wave  in  a  waveguide.   The  results  showed  a  resonance 
if  the  v/avelength  of  the  surface  wave  was  one-half  of  the 
acoustic  wavelength.   In  an  experiment  realized  in  the  Naval 
Postgraduate  School  he  verified  the  existence  of  that 
resonance. 

R.F.  Salant  [8],  in  1972,  stated  that  the  research  and 
study  of  the  theory  of  sound  propagation  in  waveguides  could 
be  and  had  been  applied  to  several  practical  problems  and 
gave  as  examples:   transmission  of  noise  through  ventilating 
ducts,  jet-engine  inlets  and  other  ducting  systems.   He 
considers  propagation  in  a  waveguide  having  two  sinusoidal 
walls  having  a  difference  of  phase,  as  opposed  by  earlier 
work  by  J.C.  Samuels  [5],  whose  waveguide  had  the  sinusoidal 
walls  in  phase. 
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Salant  concluded  that  the  phase  difference  between  the 
walls,  the  wall  wave  number  and  the  mode  number  in  the  wave- 
guide affect  strongly  the  disturbance  generated  by  the  walls. 

He  found  that  the  disturbances  when  odd  modes  are  propa- 
gating in  the  waveguide  are  extremely  different  from  the 
ones  formed  when  even  modes  are  propagating. 

A.  H.  Nayfek  [9]  in  1974  treated  again  the  acoustic 
propagation  in  ducts  with  sinusoidal  walls  using  the  method 
of  multiple  scales,  in  contrast  with  previous  works  v/here 
the  small  perturbations  method  was  used.   For  travelling 
waves,  resonance  occurs  wherever  the  v/all  wave  number  is 
equal  to  the  difference  of  the  wave  number  of  any  two  duct 
acoustic  modes. 

N.  E.  Davis  and  M.  W.  Kenyon  [1],  in  December  1974 
studied  the  effects  of  travelling  sinusoidal  surface  waves 
in  a  waveguide,  verifying  that  resonances  occur  under  certain 
conditions  relating  wavelengths  of  surface  and  acoustic  waves. 
A  waveguide  with  pressure  release  walls  was  built  and  the 
presence  of  such  resonances  verified. 

In  this  work,  a  simple  theory  is  presented  for  solutions 
of  the  acoustic  velocity  potential  in  irregular  ducts  with 
pressure  release  walls  using  the  Fourier  transform  of  the 
functions  describing  the  irregularity  of  the  wall. 
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II.   PROBLEM  DESCRIPTION 

A.   GENERAL  FORMALISM 

Consider  the  existence  of  a  waveduct,  with  transverse 

dimension  i    .   Let  the  walls  be  pressure  release  and  smooth, 
z  ^  ' 

except  for  the  upper  wall  that  will  be  irregular.   This 
waveduct  will  be  excited  by  a  monof requency  sound  source, 
and  is  assumed  to  be  of  infinite  length  in  the  x  direction. 

For  the  case  of  perfectly  smooth  pressure  release  walls, 
a  solution  for  the  wave  equation  (velocity  potential)  is, 
taking  for  simplicity  an  amplitude  of  one, 


<|)  =  sin(k   z)  cos(cot  -  k   x)  .  (2.1) 


This  represents  a  wave  travelling  in  the  positive  x  direction 
In  this  equation,  k   and  k   must  obey  the  relations: 


k   =  [(H)2  _  k  2]V2 
X      c       z 


(2.2) 


^z  ~  IT"  '   n  =  1,2 ,  . 
z 


Assume  that  the  upper  wall  of  the  waveduct  is  irregular 
and  located  at 


z   =   i^   +    e    f(x)  (2.3) 
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where  f(x)  is  the  function  that  describes  the  disturbance. 
Make  e  <<  1,  so  that  the  amplitude  of  the  disturbance  will 
be  small  compared  with  the  dimension  £  of  the  waveduct. 
Assume  a  solution  for  this  case 


^   =    ^Q    +    e    ^^   +    e"^    ^2    +    .  .  .    +    e^    (^^   +    .  .  ,  (2.4) 


In  the  above  expression,  ^^    is  as  given  in  equation  (2.1), 
representing  the  solution  for  the  undisturbed  waveduct. 

Since  the  walls  are  pressure  release,  the  boundary 
conditions  are: 


^  =  0    at 


$   =  0   at 


z  =  0 

z  =  £   +  e  f  (x) 
z 

(2.5) 
z  =  0 

^  =  ^z 


Writing  <!'  as  a  Taylor  series  expansion  about  the  point 


z  =  £^: 
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[l^   +  e  £(K)]  =  i^U^)    +  t[*i  H-  f(x)  ^]^ 


2  ^*1    1   2     ^  *0 

9z     z 


+  .  .  . 


j=0  ^• 


8^<l> 


n-3 


Sz-*     z 


+  .  .  . 


(2.6) 


With  the  imposition  of  boundary  conditions,  Eq.  (2.5) 
then  yields  the  following  set  of  equations: 


8$ 


0 


9z 


f(x) 


9$. 


8z 


f(x)  - 


1  ^'^0 


2  '    2 


I         •  f^(x) 
z 


n 


n   ^j  -  .   9^$ 
Y   f  (x)      n-3 

j=l   ^-      9z^ 


(2.7) 


In  the  assumed  solution,  $  and  ^   must  satisfy  the 

n  ^ 

homogeneous,  linear  wave  equations 

2 


r.2^       ,„2        1    9'    V 

o  $  =  (V ^  — y) 

c^    9t'' 

c   9t 


n 


=   0 


=   0 


(2.8) 
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Substitution  of  (2.1)  into  the  first  equation  of  (2.7)  gives 


=      -   k      cos  (k      z)    cos(oot   -   k      x) 


.    f(x) 


=      -    k      cos(mT)    cos(a3t   -   k      x)    f(x)     .  (2.9) 

z  X 


Equation  (2.9)  may  be  generalized  by  making  it  complex. 


h 


£ /  \  1      /   ^   i ( wt-k  x) 
=   -  f(x)  k   cos  (nir)  e      x 


(2.10) 


such  that 


Re{$ 


-i-1 


}   =   $  I 
£  1 


(2.11) 


The  solution  will  then  become 


$^  =  Re{_$j_}   . 


(2.12) 


Define : 


F(x)  =  f(x)  e"^'^^^ 


(2.13) 


The  insertion  of  definition  (2.13)  into  Eq .  (2.10)  produces 


=   -e     k   cos  (mr)  F(x) 


(2.14) 
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By  use  of  the  definition  of  the  Fourier  transform  of 
F(x) 


+00 

G(s)   =   /    F(x)  e  ^^^  dx 


and   its    inverse  (2.15) 


+00 


F(x)    =  -^      /        G(s)    e^""^    dx    , 


the  boundary  condition,  Eq.  (2 . 14),  becomes : 


-1 


ixs 
=  -e  ~  ^ 


k   cos(niT)  -^ —  f    G(s)  e     ds  .       (2.16) 

.  Z  271  ■'  — 

i  -00 


Use  as  a  trial  solution: 


00 


ia3t   ,  ,       ,      1      /-         sin(Kz)  „,    .       ixs 

(2.17) 


$T      =      -e  k      cos  (nil)    — -   J  .     ^ , . „    s       G(s)    e  ds    . 

— 1  z  217    -*         sin(K£    )       — 

_oo  2 


Since  the  above  equation  must  satisfy  Eq.  (2.8)  [the 
homogeneous,  linear  wave  equation],  the  quantity  K  must  obey 
the  relationship 


K   =   [(H)2  _  s^]^/^    .  (2.18) 
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B.   SINUSOIDAL  UPPER  WALL 

Let  the  upper  wall  of  the  waveduct  be  perturbed  in  the 
form  of  a  cosinusoide  of  amplitude  el        where   e  <<  1  . 
(See  Figs.  1  and  2)   The  upper  boundary  condition  is 


^=0   at  z   =    I    {1   +    c    cos  Yx) 


Thus  f(x)  is  specified  as 


f(x)  =  i^    cos(Yx)  .  (2.19) 


Inserting  Eq .  (2.19)  into  Eq.  (2.13) 


-ik^x 

F(x)  =  a      cos(yx)  e 
—       z 


iyx    -iyx      -ik  x 

=  I      [   ^ 1_^ ]  •  e   ^    .         (2.20) 

Z  2 


The  Fourier  transform  of  (2.20)  is: 


^z   ~    -i(k^-Y)x    -i(k^+y)x   _^ 


xs 


G(x)  =  -^  /    [e    "     +  e    ^    ]  e  "'^^  ds 


2 

—00 


=  T\Z^    [6(s+k^-Y)  +  6(s+k^+Y)]  (2.21) 


and  the  solution  will  be: 


—  00  7.  <  ' 


•   e^^^  ds  .   (2.22) 
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tfx 


■"+X 


FIGURE^ 1.   The  Unperturbed  Duct 


t+1 


+x 


FIGURE. 2.   The  Perturbed  Duct 
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Using  the  properties  of  the  6-f unction,  the  integral  is 
readily  evaluated  so  that 


$,  =  -e    k    cos  nTT 

— 1  2  Z 


r  /  Wx  2     ,.       ,  X  2  , 1/2 
r      ^^^^^c^   -  ^V^^  ^     ^     i[a3t-(k^+Y)x] 
_  \ e  + 


sin[(^)^  -  (k  +Y)2]l/2  ^ 


r  /  Wv  2    ,,     V  2  ,  1/2 
sin[(-)   -  (k^-Y)  ]     z     i[^t-(k  -Y)x]  ^ 

+   e  y    .  (2.23) 

r  /  tOx  2     ,,     X  2  T  1/2  0 
sin[(-)   -  (k^-Y)  ]     ^2 


Since  <|),  =  Re  {_$,  }  ,   the  first  order  correction  i; 


$^  =  -k  £   cos  (nrr) 


.+ 


1  f  sin  K  z      r    .  /T  .  \  1  . 

2  i  sin  K-£   ^°^f"^  -  (V^^^^  -^ 


^   sin  K-£   cosl^t-  (Vt)x]}    ,  (2.24) 

z  ^ 


where   K"  =  [  (^)  ^  -  (k^±Y)^]  "^^^       .  (2.25) 


1.   Validity  and  Confirmation  of  Solution 

The  above  solution  is  consistent  with  the  results 
presented  by  N.  E.  Davis  and  M.  W.  Kenyon  [1].   The  waveguide, 
in  that  case  is  perturbed  by  a  traveling  surface  wave  of 
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amplitude  A.   The  upper  boundary  condition  is  $  =  0   at 
z  =  il  +    e   F(x,t)  ,   where 


F(x,t)  =  -£   cos(fit-yx)    and    e  =  -j — 

z 


By  making  Q,   =    0      in  Eq.  (2.12)  of  [1]  and  discarding  the 
functional  dependence  in  the  y-direction,  the  result  agrees 
completely  with  Eq.  (2.22)  . 

2 .   Properties  of  the  Solution 

In  Eq.  (2.24)  ,  the  propagation  constant  in  the 
z-direction  is 


K-  =  [(^)^  -  (k  ±y)^]-^^^   .  (2.26) 

O  X 


/-I  •      / (jO\  2    ,  2  ,  ,  2 
Since   (— )   =  k   +  k    , 
c      X      z      ' 


k 

(1  ±  2  -^)   ,/^ 

K-  =  k   [1 :^   \^^^       .  (2.27) 

^        (k/Y)'' 


If   K    is  real  and: 


K*  £    =   miT    ,    m=l,2,...,  (2.28) 


the  conditions  for  a  resonant  situation  exist.   From 

nTT 
Eq.  (2.2),   k   =  - —   ,   n=l,2,...   .   Inserting  this  into 

^    ^z 
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Eq.  (2.28)  the  following  expression  results,  which  states 
the  relation  necessary  for  a  resonance: 


or 


(1  i  2  yv) 

n[l  -  5 ]  '         =     m 


,m,2      ,    (1*2  Vy) 


O    =   1 ? .  (2.29) 


The  resonance  curves  are  shown  in  Figure  3  as  a  plot  of   k  /y 

2 
versus   (k  /y)   .   This  is  in  accordance  with  Eq.  (2.16)  and 


Figure  2  of  [1] ,  and  [3] . 
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^ 

^-^y 


FIGURE  3.   The  Re 


sonance  Curves 
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III.   ANALYSIS 

Using  spectral  analysis  it  is  possible  to  find  the 
solution  to  the  problem  involving  any  special  boundary 
condition  related  to  the  upper  wall  of  the  waveduct .   The 
only  thing  that  has  to  be  known  is  the  Fourier  transform  of 
the  combination 


-ik  X 
f(x)  •  e    ^      ,  (3.1) 


v/here  f  (x)  represents  the  function  describing  the  upper  wall 
of  the  waveduct. 

Knowing  this,  the  problem  is  reduced  to  solve  the 
integral 

_    1    f      sm  Kz    ^,    .       ixs  ,  ,-,  ^. 

I  =  TT—      — = TTT—   G(s)e     ds,  (3.2) 

27T    J       sm  KZ        —  ' 

-oo  2 

where  G(s)  is  the  Fourier  transform  of  (3.1) , 


°°         -1  (k  +s)  X   ,  /  o  o\ 

r,  I    \  (      HI    \  X       dx    .  (3.3) 

G(s)  =   J   f(x)  e 


In  the  problem  treated  and  solved  before,  the  upper  wall  was 
described  by  a  very  simple  function  and  the  evaluation  of 
the  integral  did  not  pose  any  special  difficulties.   But,  if 
f(x)  is  a  more  complicated  function,  a  careful  analysis  of 
(3.2)  is  necessary  in  order  to  solve  it. 
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A.   THE  SOLUTION  INTEGRAL 

If  s  is  generalized  to  a  complex  quantity  (s  is  the 
variable  of  integration) ,  a  careful  study  of  the  behavior  of 
the  function  under  the  integral  sign  in  (3.2)  in  the  complex 
s-plane  is  necessary. 

The  function  to  study  is: 


sin[(-)   -  s  ]  '^   z 
f(s)   =   G(s)  e^^^   .    (3.4) 


r  ,  COv  2        2  T  1/2  n 

sin[(— )   -  s  ]  ^   £ 
c  z 


1.   The  Poles 

Assuming,  for  the  moment,  that  G(s)  does  not  have 
any  poles,  interest  is  focused  in  the  ratio  of  sines  in  (3.4) 
Since  [2] 


2 

sin  z  =  z   n    (1 ~j  )         ,  (3.5) 

m=l       m  77 


f(s)  can  be  written  as 

[(^)2-s2]  z^, 

ooi  C  ) 

r,a3v2      2,1/2              „       ),  {„,    x^ixs 

[(-)    -s    ]    /       z         n        1 ^-^ G(s)e 

f(s)    =   JB^^-J fil- 1 (3.6) 

[(^)2-s2]l/2    ,^      I      i  [(§)^-s^]    C    I 

^  ^   m=l    )  ^2    2  ( 

^  m  77  ' 


or   simply 
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CO  [11  r/C0\2  2,  2 

f(c\    -   Jl.         n         -^m    TT      -     [(-)     -s    ]     z     ,  . 


f    2    2 

|m    IT      - 


c  ■■     ^z 


There  are  no  branch  cuts,  the  poles  of  f (s)  are  simple  and 
located  at 


s^  =  ±  [(^)^  -  (j^)^]  ^^^  (3.8) 

2 


m= 1,2,... 

It  is  possible  to  arrive  at  the  same  result  by  simply  looking 
at  the  denominator  of  (3.4).   The  poles  will  be  located  at 

j^^j2_g2jl/2  ^   ^  ^  ^^  (3^9) 

m=l, 2 , . . . 

Manipulation  reveals  that  Eq.  (3.9)  is  equivalent  to 
Eq.  (3.8),  but  here  it  is  not  evident  that  the  poles  are 
simple,  and  that  there  are  no  branch  cuts. 
From  Eqs .  (2.2)  it  is  known  that 


(iii)2  =  k  2  ^    2 
C        X       z 


and 
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k   =111   .   n=l  2 


Using  this  in  Eq.  (3.8) ,  it  can  be  shown  that  the  location 
of  the  poles  are  at 


^m  "  -  f^x^  "^  (■^)^(n^-in^)]^^^  (3.10) 

z 


n , m= 1,2,... 

Looking  at  this  expression,  it  can  bee  seen  that  the  poles 
will  be  located  in  the  real  axis  if 


k^^  >  {~)^{n^-m^)  (3.11) 

z 


n , m= 1,2,...   , 

and  on  the  imaginary  axis  if  the  inequality  is  reversed. 
Let  us  choose  n=l  for  all  that  follows.   This  means 
that  only  the  lov/est  mode  can  propagate  energy.   It  can  be 
seen  that  as  m  increases  the  associated  poles  are  found 
successively  closer  to  the  origin  on  the  real  axis  until 
above  some  critical  value,  m  ,  the  poles  lie  on  the  imaginary 
axis  with  distance  from  the  origin  increasing  with  m.   The 
first  pole,  s, ,  will  be  located  at 

Sn  =  ±  k    .  (3.12) 

1      x 
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If  it  is  assumed  that  m  =  3  ,   the  fourth  pole  is  already 
on  the  imaginary  axis.   In  Table  I  and  Figure  4  a  situation 
like  this  is  shown. 

2 .   Effects  of  Absorption 

If  the  effects  of  absorption  are  introduced,  the 
wave  equation  that  governs  the  propagation  inside  the  duct 
is  modified  and  becomes: 


(°'  -  2  ^It'  *n   =   °  '3.13) 


n  =  0,1/2,...   , 


where  the  a  's  are  the  coefficients  of  absorption  responsible 

for  losses.   The  a  's  can  be  assumed  as  having  the  same 

n  ^ 

magnitude  since  all  the  $'s  have  identical  form,  and  the 
propagation  constant  in  the  x-direction  will  become  complex. 
It  is  known  that  the  $'s  have  the  form 


i [wt-k  x] 
sin(k^z)e      ^  (3.14) 


Using  this  in  Eq.  (3.13) ,  the  value  of  k   is 


k   =  k   -  ia  -~-   =  k   -  i3 
— X    X       ck     X    ^ 

X 


where 


(3.15) 


3  =  a   '^ 


ck 

X 
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TABLE  I 
The  Poles  of   f(s) 


n=l  ,  m=l   :  St  =  ±  k 

'  1       X 


2    -,  TT  .  2,1/2 


n=l  ,  in=2   :  ^2  "^  -  ^^x   "  ^^f"^^ 


n=l  ,  m=3   :  ^3  =  1  [k^^  -  Si^-)^]^^^ 


n=l  ,  m=4   :  5^  =  1  [k^^  -  15  (-^L)  2  j  1/2 

z 


=   ±  i  [15(^)2  -  k^2,l/2 
z 


n=l  ,  m=5   :  g^  =   ±  i  [24  (■^)  ^  -  k^^]^^^ 

z 


n=l  ,    m=6,7,...  :     s^  =   ±  i  [  (m^-l)  C-^)  ^  -  k^^l^^^ 

z 
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FIGURE    4.      The   Poles    of    f 


(s) 
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is  the  total  absorption  coefficient  in  the  x-direction  and 
will  be  considered  very  small. 

It  is  necessary  to  know  what  is  the  effect  of  the 
introduction  of  absorption  in  the  location  of  the  poles  of 
f  (s)  ,  especially  the  ones  located  in  the  real  axis.   If  s-, 
is  located  at   s,  =   ±k   ,   it  can  be  seen  that  the  poles 
in  the  negative  part  of  the  real  axis  will  be  moved  up  and 
the  poles  located  in  the  positive  part  will  move  down 
(Figure  5)  .   We  may  nov;  choose  a  contour  in  order  to  solve 
the  integral. 

B.   THE  CONTOUR  INTEGR?^L 

The  integral  to  be  solved  is 

00 

I   =   /   f(s)  ds  (3.16) 

—  00 

where  f(s)  is  the  expression  in  (3.4). 
1.   Choice  of  the  Contour 

The  integral  in  (3.16)  is  very  similar  to  the  one 
Isakovich  used  in  his  work  [4].   It  reduces  to   2iTi   times 
the  summation  of  the  residues  at  the  poles  of  f (s) . 

The  contour  to  consider  will  consist  of  the  segment 
[-R,+R]  on  the  real  axis  and  the  upper  semicircle  C_,,  as 
shown  in  Figure  5. 

Therefore 


«o  +R 

I  =  /  f  (s)  ds  =  lim   /     f(s)  ds  =  2i\i    J  Res^^  f(s)  .    (3.17) 

-oo  R-^oo  -R  m 
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FIGURE    5. 


or?.^e^%o^^ef  ^r?lS-   -   -   -auon 
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simplify   Eq.     (3.4)    by   defining   H(s)    as: 


r  /  COv  2  2  ,  1/2 

sin[(— )       -    s    ]    ^         z 

H(s)    = .  (3.18) 

r    /  OOv    2  2  ,    . 

sin[  (-)       -    s    ]  I 


For  the  case  where  G(s)  does  not  have  any  poles  the 
integral  will  be: 


00  ixs 

I  =  /  f(s)  ds  =  2TTi   y   G(s  )  e   ^  Res{H(s  )}    (3.19) 
-00  m 


On  the  other  hand,  if  the  Fourier  transform  of  F(x)  has  poles 
of  its  own,  at  s  =  s  ,  where 

g 


-^  0   ,  (3.20) 


£%) 


the  expression  to  evaluate  the  integral  is  modified  as 
follows : 


00  ixs 

I  =  /   f(s)  ds  =  2Tri  )      G(s  )  e   ^   Res{H(s^)} 
;  ^   —  ^  m  m 

-00  ni 

ixs 
+    27Ti      y      H(s    )     e         ^   Res{G(s    )}  (3.21) 

g  ^ 


In  either  case,  the  poles  of  interest  will  be  the  ones 
included  within  the  contour  chosen. 
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2.   Definition  of  the  Transform  of  F(x) 

In  order  always  to  be  able  to  close  the  contour  in 
the  upper  half  plane,  the  definition  of  the  Fourier  transform 
of  F(x)  must  be  altered. 
We  define 


G(s)  =  /  F(x)  e'"-''"''  dx 


with  inverse 


CO 

F(x)  =  ^  /  G(s)  e-^""^  ds 


27r 


The  upper  sign  of  the  exponent  will  represent  the 
positive  x-part  of  the  waveduct  with  waves  standing  or 
propagating  to  the  right  and  the  lower  sign  will  represent 
the  negative  x-part  of  the  waveduct  with  waves  standing  or 
propagating  to  the  left. 
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IV.  SPECIAL  CASES 

In  order  to  simplify  the  problem  and  apply  this  theory 
to  some  specific  cases  we  shall  impose  some  conditions  on 
the  problem.   Since 


(M)2  =  3^  2  ^  j^  2  ^   ^ 

C        X       z     ' 


in  the  duct,  the  cut-off  frequency  will  be 


cj^  =  c  k^   .  (4.2) 


The  sound  source  that  is  going  to  excite  the  duct  will  have 
a  monofrequency  of 


03  =  A  c  k    ,  (4.3) 

z 


where  1  <  A  <  2  (4.4) 

to  assure  propagation  only  in  the  lowest  modes.   Under  this 
condition  and  knowing  that 


K  =  Tl    ■  <^-5> 


(only  lo^^estmode  allowed  in  z,    n=l)  the  location  of  the 
poles  of  f(s)  will  be  at 
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It  can  be  seen  that  only  one  pole  will  be  on  the  real 

axis  (s-|)  ,  and  that  all  others  are  already  on  the  imaginary 

axis. 

The  general  expression  for  the  location  of  the 
poles  of  f (s) ,  assuming  G(s)  has  no  poles  is: 


s^  =  ±[kj   +    (^)^(l  -  m^)]^'^^  (4.6) 

m        X      x,^ 


and   the  m        residue   is 


.     ,  r  ,a3v  2      2,       2    \    '^  ^ ,    V       ±ixs 

sin-{  [  (g)    -s    ]    z       }  •    G(s)    e 


cosi[(^)2-s2]    .^2    JV2    ^    1    |^^^^2_^23    ^z'P'-^^^z' 

(4.7) 


evaluated  at  s  =  s_  . 

m 


This  will  become,  after  manipulation, 


_       /  T .  m+l     z   .   miT    _.  /   V   ±ixs^  ,.    o\ 

m  z      z 


m  =  1,2,... 
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The   effect  on  acoustic  velocity  potential  due  to  the  m 
pole  is  [using  Eq.  (2.17)] 


th 


itot 


1  r 

k  TT-  *  2iTi   )   Re5 


'z  211 


m 


m 


(4.9) 


If  G(s)  has  poles  of  its  own,  (4.9)  will  be  modified 
and  the  effect  on  the  velocity  acoustic  potential  will  be 


loot  ,    1  .  ^  . 

Z  27T 


I      Res 


m 


m 


+  I      H(s  )  e 

g    ^ 


±ixs 


^  Res  ] {G(s  )  } 

gl   g 


(4.10) 


A.   THE  6-FUNCTION  DISTURBANCE 

Define  f(x),  the  function  describing  the  wall  disturbance 
as 


f(x)  =  -Z^    6(x) 


(4.11) 


Then 


F(x)  =  -£   5  (x)  e 
—        z 


-ik  X 

X 


(4.12) 


and 


G(s)  =  -I 


(4.13) 


with  no  poles. 
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The  poles  of  f(s)  will  be  located  according  to  Eq.  (4.6) 
The  first  pole,  the  only  one  in  the  real  axis  will  be  at 


s^  =  -k^  (4.14) 


and  the  second  one,  on  the  imaginary  axis  will  be  at 


S2  =  +  i  [(^)^(22-l)  -  k^^]^'^^        (4.15) 


=  +  iP 


where  P  is  defined  as 


P  =  [3(^)^  -  k^^]^/2   .  (4.16) 

z 


The  overall  situation  is  shown  in  Figures  6  and  7.   The 
residue  of  f(s)  at  s,  ,  Res-|(s-,)  will  be,  using  Eq.  (4.8), 

^z  +^^x^ 

Res*  (s,)  =  ?— -  sin  k  z  e  (4.17) 

X 

The  acoustic  potential  due  to  this  residue  will  be  calculated 
using  (4.9) ,  and 

2 

r  "v      k  i(a)t;k  x) 

^^|ReSj^(Sj^)  y   =  "Y"  sin  k^z  i  e      ^     .       (4.18) 
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t^ 


FIGURE  6.   The  6-Function  Disturbance 


FIGURE  7.   The  Poles  of  Interest 
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since  $,  =  Re{  $_,  }  ,      in  this  case 


k^2 

<|)n  =  -  -y^ —  sin  k^z  sin((jot+k^x)  (4.19) 


The  above  equation  represents  two  scattered  traveling  waves 
emanating  from  x=0,  one  going  in  the  positive  x-direction  and 
the  other  in  the  negative  x-direction. 

The  residue  of  f(s)  at  S2 /  ReSp(s„),  will  be,  using 
Eq.  (4.8), 


2k 

Res2(s2)  =  -i  -^   sin  2k  z  e"^  (4.20) 


with  a  corresponding  acoustic  velocity  potential  of 

2k  ^  •  4-   -  p 

$^[Res2(s2)]  =  — ^  sin  2k^z  e^^"^  e^^    .      (4.21) 

This  represents  two  disturbances  in  the  acoustic  field, 

decaying  in  the  ± x-direction.   It  can  be  seen  that  far  away 

from  the  wall  disturbance  their  effect  will  be  negligible. 

Very  close  to  the  wall  disturbance,  x=0,  the  effects  will 

be  considerable  and  depend  on  the  size  of  P_/  diminishing 

^  m 

as  P_  increases  with  m.   P_  is  defined  as 
m  m 


^m  =  [(m^-l)(-^)^  -  k^2]^/2   ,    m=3,4,... 
z 


However,  even  at  great  distances  from  the  irregularities 
in  the  wall  the  acoustic  field  is  affected  by  $,  [Res-,  (s, )  ]  . 
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The  total  field,  discarding  the  effects  of  decaying 
disturbances  as  in  (4.21)  is: 


k  2 

=  sin  k  z  [cos  (tot-k  x)  -  e  -r—-   sin(cot+k  x)  ]       (4.22) 
z  ^        ^x  ^ 


Note  that  the  traveling  perturbation  is  9  0**  out  of  phase  with 
the  initial  wave. 

B.   THE  DECAYING  WALL 

Define  f(x),  the  function  describing  the  wall  disturbance 
as 


f  (x)  =  -l^   e"i^^l  (4.23) 


Then 


and 


_ ,   ,   -ik  X 
F(x)  =  -£^  e  l^^l  e    ^  (4.24) 


2^z 
G(s)  = (4.25) 


|a|[l  +  (|  +  k^)2] 


a    X 


with  poles  located  at 


s„  =  a(-k^  ±  i)  (4.26) 
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The  poles  of  interest  of  f(s)  will  be  located  at 


s,  =  -k 
1      X 


S2  =  +iP  (t  .27) 


s   =  a(-k  +  i) 
g        X    ' 


The  overall  situation  is  shown  in  Figures  8  and  9.   Since 
G(s)  has  poles,  Eq .  (3.21)  must  be  used  to  evaluate  the 
integral. 

The  residue  of  f(s)  at  s,,  Res, (s,)  will  be 


2   k  2  ik  X 

^^^l^^l^  =  a  ir  ^^^  ^z^      2  ,  ^%,   ,,2  ^    ""       (^-28) 
X         a  +  k   (a-1) 

X 

with  a  corresponding  acoustic  velocity  potential  of: 

k^^       a  i(wt+k^x) 

_$-|_  [Resj^  (Sj^)  ]  =  i  2  -j^ 2 2 2  ^^^  ^z^  ® 

X  a  +  k   (a-1) 

^  (4.29) 


representing  two  traveling  waves  emanating  from  x=0,  going 
in  opposite  directions,  and 

2k  ^ 

^^   =   Re{^^}   =  -  ^1^ 1 2"  sin  k  z  sin(wt+k  x) 

X  a  +  k   (a-1) 

^  (4.30) 
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I 


FIGURE  8.   The  Decaying  Wall 


FIGURE  9.   The  Poles  of  Interest 
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The  residue  of  f(s)  at  S2/  Res2(s2)  will  be: 

z  +Px 

Res^(s^)  =  r— ^—  sin  2k  z  e  (4.31) 

iPa[l  +  (^  +  k^)^]        ^ 
with  a  corresponding  acoustic  velocity  potential  of 

4k  2 

.    rr^      /     \  1  Z  •     m         itOt    +PX 

^-t   [ResT(So)  J  =  r— :r—   Sin  2k  z  e     e 

-'     2   2      Pa[l  +  (i^)  +  k  )2] 

^      ^  (4.32) 


This  represents  two  standing  disturbances  decaying  with 
X  and  negligible  at  great  distances  from  the  origin. 

The  residue  of  f (s)  at  s   will  be  computed  using  (3.21) 

y 


and  is 


±ixs 
H(s  )  e     ^  Res{G(s  ) }  (4.33) 

y  y 


or 


sin  v^TIc   z_  ^^ax  ^^fV  .  ^  ^^^^^^ 


sin  /B+iC   £ 
z 


z 


where: 


2      2       2      2 

B  =  k  "^  +  k  ^(1  -  a  )  +  a^ 
z      X 


C  =  2a^k 

X 


(4.35) 
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The  corresponding  acoustic  velocity  potential  is: 


sin  /B+iC   z    _    i(wt+k  x) 

$,[Res(s)]=k    e^^""  e      ""  (4.36) 

^  ^      ^  ^  ^  sin  /B+IC    I 

z 


This  corresponds  to  two  traveling-decaying  waves  (both 
in  X  and  z  directions)  going  in  opposite  directions  from 
x=0 .   At  large  distances  from  the  origin  its  effect  will  be 
practically  nil. 

Combining  results,  only  valid  at   x  >>  0   or   x  <<  0  , 
we  have: 


<I>    =    'J>Q    +    G    (I', 


2k    ^ 

=   sin   k    z    [cos(ajt-k   x)    -    e   —r- -^ ^ s-  sin(wt+k    )] 

X      a    +  k      (a-1) 

X 

(4.37) 


It  can  be  noticed,  again,  that  the  traveling  disturbance 
is  90'*  out  of  phase  with  the  original  propagating  wave. 
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V.   RESULTS  AND  CONCLUSIONS 

The  results  obtained  show  that  the  effects  of  very  small 
irregularities  in  the  physical  boundaries  of  the  duct  can 
have  appreciable  importance  in  the  waveduct  sound  propagation 
even  at  very  large  distances  from  the  disturbances.   In  the 
case  of  the  sinusoidal  wall,  resonances  can  exist,  which 
could  have  the  same  order  of  magnitude  of  the  original 
propagating  wave.   These  resonances  do  not  appear  in  the 
two  other  special  cases  studied,  the  6-function  disturbance 
and  the  decaying  wall.   It  seems  that  the  theory  works  and 
can  be  applied  to  other  types  of  corrugations,  once  the 
Fourier  transformation  of  the  function  describing  the  wall 
disturbance  is  known.   However,  complications  can  arise 
depending  on  the  complexity  of  the  above  mentioned  Fourier 
transform,  making  the  solution  integral  difficult  to  evaluate, 
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